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Abstract 

The general relativistic accretion onto a black hole is investigated in 
which the motion is steady and spherically symmetrical, the gas being at 
rest at infinity. Two models with different equations of state are compared. 
Numerical calculations show that the predictions of the models are similar 
in most aspects. In the ultrarelativistic regime the allowed band of the 
asymptotic speed of sound and the mass accretion rate can be markedly 
different. 



1 Introduction 

The accretion of gas on compact objects (white dwarfs, neutron stars, black 
holes) has not been entirely investigated by now, even in the simplest case of 
spherically symmetrical systems. The study of accretion has its beginnings in 
the paper presented by Bondi (1952) pQ. He considered spherically symetrical 
accretion on the basis of Newtonian gravity. Further progress has been made 
by Michel (1972) [2] and Shapiro and Teukolsky (1983) [Hj who gave a general 
relativistic version of the Bondi model (the (p — n) model). Another relativistic 
generalization was given by Malec (1999) gj (the (p - p) model). 
It is not clear which equation of state is appropriate in the description of rel- 
ativistic collapsing gas. There are two commonly used polytropic equations of 
state: p = Kp T [3j and p = Cn T Here p is the pressure, p is the density 
and n is the baryonic mass density. The intention of this paper is to compare 
predictions of both models concerning the sound velocity, fluid velocity, density 
and mass accretion rate. Here the {p — p) model and the (p — n) model denotes 
the model with the equation of state given by p = Kp T and p = Cn T , respec- 
tively. 

We also show that the {jp — n) model gives an upper bound for the asymptotic 
speed of sound [H]- 

The order of this work is as follows. In Sections 2 and 3 we briefly present 



1 



(p — n) and (p — p) models. Section 4 is dedicated to the derivation of the after- 
mentioned limit on a 2 ^ and a 2 s for the (p — n) model. In Section 5 we compare 
predictions of both models using the results of numerical calculations. 
In the course of the paper we set G = c = 1 everywhere. 



2 The (p — n) Model of Stationary Accretion 

Here we shall give a short briefings of this model following Shapiro and Teukolsky 
PJ ■ A polytropic equation of state is assumed 

p = Kn T , (1) 

where K and T are constant. The velocity of sound is given by 

9 dP 

The boundary conditions are as follows: the gas at rest is described by the 
baryonic mass density n^ and total energy-mass density p^. We omit details 
(can be found in [3]) and write down the final equations. One can find that at 
the sonic point R = R s the speed of sound a and the infall velocity u satisfy 
relation 
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From the relativistic Euler equations in Schwarzschild coordinates we get 
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'-T^)('^JH 1 + f^J' (4) 

Making use of conservation of the baryonic mass: AirnuR 2 = const, and rear- 
ranging Eq. one finds that 

1 (air 2M~~X^ (l + 'ial\ V2 _ 

u= AR*{t*r--[i +u ) \-^t) ■ (5) 

Inverting both sides of Eq. and evaluating at the sonic point with the aid 
of Eq. Q one gets the key equation for our considerations: 
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Employing Eq. |jSJ one can describe mass accretion rate by [3] 

^ /l + 3a^ 3/2 



7\/ = ^M 2 n 0O (^.yTT3^) (^P) ■ (7) 
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3 The (p — p) Model of Stationary Accretion 



Here we describe the relativistic model shown in 0j. A suitable choice of an 
integral gauge condition leads to the comoving coordinates formulation that is 
particularly suitable for the description of self-gravitating fluid. Spherically 
symmetric line element is given by 

ds 2 = -N 2 dt 2 + adr 2 + R 2 dQ 2 + R 2 sin 2 &d(j) 2 , (8) 

where iV, a and R depend on t (asymptotic time variable) and the radius r. 
The energy-momentum tensor of self-gravitating fluid in comoving coordinates 
is given by 

Tp V = (p + yju^Uy + pg^ , (9) 
where = -1. Notice that p = T r r = Tf. 

The rate of mass accretion M along orbits of a constant areal radius R is equal 

to m 

M(R) = -4ttNR 2 u( P + P), (10) 

where 

u = d R/N. (11) 

The effect of backreaction is neglected, that is the change of geometry caused 
by infalling gas is regarded to be negligible. Then one finds 0] that 



W .f = ^l-f + „*, (12) 



where p is the mean curvature, and 



2R S ) 3 2R S v ; 

must hold in the sonic point. Assuming that the equation of state is given by 

P = I<P F , (14) 

where the constant T belongs to the interval 1 < T < | , and defining (as usual) 
the velocity of sound as a 2 = d p p one arrives at 0] 

where k = (T — l)/T and the integration constant a 2 ^ is equal to the asymptotic 
velocity of sound at infinity. 
One can find that 



u 2 
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1 + 72 • (16) 
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It should be emphasized that Eqs. lj*H?|l . (|T5|) and (|Po|) form a purely alge- 
braic system of equations describing the fluid accretion in a fixed space-time 
(Schwarzschild) geometry. 

From the relation between pressure and energy density one obtains that 



P = P 00 (a/a oa ) 2/(r 1} = Poo ( — n- + 



1 



i/(r-i) 



(17) 



where the constant poo is the asymptotic mass density of the collapsing fluid. 
Substituting Eqs. (jl3|) and 1)15(1 into rearranged l|10|l we find that the mass 
accretion rate can be described by means of the formula |1] : 



(5-3I-) 

2(r-i) 



(18) 



4 Limit on in the (p — n) Model 

It will be useful to transform Eq. into 



l + 3ai) = 



i-<&/(r-i) 



-i 2 



i-a2/(r-i) 



(19) 



Let us assume that < a 2 s < a 2 lax and X = (1 + 3a^ aa ,). Thus values of the 
left hand side of Eq. i|19|) belong to the range [1;X]. Hence 



1 < 



Ll-a2/(T-l) 



< X. 



Solving (J20J) we get 



< a 2 s < (T - 1) 1 



■ CXI- 

X 



and 



2(T-l)-a 2 0O >a 2 s >(T-l) 1 



(20) 



(21) 



(22) 



These conditions can give us a restriction on the asymptotic velocity a 2 ^ . Simple 
calculations lead to the inequality [5] 



«L< (r-i). 



(23) 



5 Numerical Calculations 

In this section we compare both models numerically referring to certain param- 
eters important for the description of the process of accretion. 
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5.1 Evaluation of Parameter a 2 s 

First we analyse formula ^(a 2 ^) in the (p— ri) model given by Eq. ©. The ana- 
lytical bound on a 2 ^ shown in previous section has been confirmed by numerical 
calculations. In Fig. ^ we show all solutions of Eq. © for Y = 1.3. As one can 
see there are four families of solutions corresponding to various combinations of 
the branches (1) - (4). We eliminate from our considerations the branches (3) 
and (4) which describe the case a 2 < a 2 ^. It is not clear that the branch (2) 
should be rejected. However it seems peculiar for us that the asymptotic speed 
of sound decreases while the speed of sound in the sonic point increases (e.g., 
for a 2 ^ =0, a 2 s has a nonzero value). Notice that this branch is not represented 
at the Newtonian level since in the Bondi model a^(a^ = 0) = 0. In the (p — n) 
model it is assumed that for at large R the flow satisfies condition u <C 1 and 

1 i 1 1 1 1 1 
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Figure 1: Plot of a 2 in terms of a 2 ^. All solutions are shown. The branches (3) 
and (4) have no physical meaning (a 2 < a 2 ^). In this paper only the branch (1) 
is being analysed. 

is subsonic with u 2 < a 2 (e.g. as R — ► oo, u — » 0, a — * Ooo). For the branch (2) 
numerical calculations of the behaviour of parameter u as a function of R yield 
relativistic values of the parameter u, even if R is very big (Tab. 2J. Therefore 
we exclude the branch (2) from what follows. Next we analyse the solutions 
of Eq. |JBJ for certain values of T and compare them to the solutions of Eq. 
<|15[) in the (p — p) model. Our calculations are shown in Fig. [21 Numerical 
results confirm previous analytical estimations of the cut-off of the parameter 
a 2 . As one can see in the (p — n) model the greater Y is the greater value of a 2 
is reached. In contrast no such limitation appears in the (p — p) model. 
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Figure 2: Plot of a 2 in terms of for three different values of the adiabatic 
index: T = 1.1. T = 1.2. and T = 1.6. Dotted and solid curves refers to the 
(p — n) model and the (p — p) model, respectively. 



u(R) [c] 


R [M] 


0.981071 
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0.873213 


10 


0.776209 


50 


0.763217 


10 2 


0.751332 


10 2 


0.750134 


10 4 


0.750014 


10 5 


0.750002 
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Table 1: The values of u(R) for the supercritical branch (2). 
5.2 Fluid Velocity 

Fluid velocity as a function of a distance is described by Eqs. JSJ and 1161) for 
the (p — n) model and the (p — p) model, respectively. It rises monotonically 
as the radius tends to the event horizon. Comparing both models we assume 
the same asymptotic sound velocity a 2 ^. As one can see (Figs. |3]and0J both 
models predict similar values of u. 

We noticed that the greater T the slower fluid velocity is at the given distance 
R. 

Next conclusion is the confirmation of the fact (previously stated in 01) that 
the value of u near the horizon strongly depends on the location of the sonic 
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point R s . The further the sonic point the larger the fluid velocity and the closer 
to the speed of light at R = 2M. 
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Figure 3: Plot of fluid velocity u as a function of a radius R for T = 1.1 and 
= 0.099. Dotted and solid curves refer to the (p — n) model and the (p — p) 
model, respectively. 



5.3 Density Profile 

We recall here that the main difference between the (p — n) model and the (p — p) 
model lies in the equations of state: p = Cn T and p = Kp r , respectively. We 
can relate n and p by 

n = exp (j dp JT^) (24) 

that can be integrated with the result 



p(l + hp ) =p[l + — 



i/(r-i) 

r 

Given the (p — p) polytropic model one can always find n. And conversely, 
one can find p, given the polytropic (p — n) model (Hj- The preceding equation 
yields rioo = p^ if <C 1; the same is true in the alternative description 
(n — p — > p) under the condition Poo/(T — 1) <C Poo®- According to the 
numerical calculations when matter approaches the horizon its density increases. 
We also noticed that the location of the sonic point R s plays a very important 
role. If it is situated close to the horizon the density of matter there increases 
approximately to 10 x p M . On the other hand if R s ^> 2M the density of 
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Figure 4: Plot of fluid velocity u as a function of a radius i? for T = 1.4 and 
a^, = 0.099. Dotted and solid curves refer to the (p — n) model and the (p — p) 
model, respectively. 
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Figure 5: Plot of nondimensional density profile p/poo as a function of the 
distance R for the (p — p) model. The asymptotic velocity = 0.099 for both 
r = 1.1 (solid curve) and T = 1.4 (dotted curve). 

matter approaching the horizon becomes few orders of magnitude greater than 
the asymptotic density. The predictions of the two models agree in the full 



spectrum of the values of index T. 




Figure 6: Plot of nondimensional baryon density number profile n/rioo as a 
function of the distance R for the (p — n) model. The asymptotic velocity 
= 0.099 for both V = 1.1 (solid curve) and Y = 1.4 (dotted curve). 



5.4 Mass Accretion Rate 

In this subsection we compare the most important parameter to the description 
of accretion: mass accretion rate M. For simplicity we introduce the parameter 
Q which is defined as the ratio of mass accretion rate in relativistic model and 
the mass accretion rate predicted by the Bondi model: M = HAIb- Hence fi 
can be interpreted as relativistic correction factor. 
In the (p — n) model this parameter, with help of J7J), is expressed by 



« = <4 -Wl 



3a 



l + 3ai 



3r 



(25) 



while in the (p — p) model using (|18|l we get 



n = 



(5-3r)a f^^ 
2aL 



l + 3af 1 



(26) 



The comparison of the parameter fl for the two models (Figs. ElandEl leads 
to the conclusion that they slightly differ in a full range of allowed a^, but it 
should be emphasized that the accretion in the (p — n) model is more efficient. 
Next we compare the relativistic correction factors as functions of the adiabatic 
index. We consider here an ultrarelativistic regime, i.e. we assume the maximum 
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possible value of a 2 . 

In 0] it was shown that for the (p — p) model the relativistic correction factor 
satisfies 

4fl + ~J>n>1.6(l + ~). (27) 

We confirm here that the values of a parameter £1 belong to the range defined 
by H27|) . However we revealed earlier [Sj, that the factor is not a monotonic 
function of T and for T 1.46 it has a minimum of a value Q 4.77. This 
is again confirmed by the present calculations. For the (p — n) model Q rises 
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Figure 7: Plot of the relativistic correction factor f2 as a function of the adiabatic 
index T for both models in ultrarclativistic regime. For each T the maximum 
possible a 2 s is set. 

monotonically as T increases (Fig. 0). 

It should be mentioned that in nonrelativistic case a 2 <C 1 the relativistic correc- 
tion factor is close to 1 (Fig. 8), in full agreement with theoretical expectations 
□ 

6 Conclusions 

We examined two models of stationary and spherically symetrical accretion of 
gas onto a black hole. We show that both models essentially agree as it concerns 
quantities such as fluid velocity u, density profile and the mass accretion rate 
M. 

What drastically differs the models is the bound on the sound velocity a 2 s which 
has been found in the (p — n) model. No such restriction appears in the (p — p) 
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Figure 8: Plot of relativistic correction factor fl as a function of asymptotic 
velocity of sound for fixed value of T = 1.1. Dotted and solid curves refers to 
the (p — n) model and the {p — p) model, respectively. 
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Figure 9: Plot of the relativistic correction factor 11 as a function of asymptotic 
velocity of sound for fixed value of T = 1.4. Dotted and solid curves refer to the 
(p — n) model and the (p — p) model, respectively. 

model. It makes this model more advantageous especially for the values of 
adiabatic index T close to 1 where the (p — n) model provides the solutions 
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only in a very narrow range. Another interesting difference can be observed in 
the ultrarelativistic regime and for T close to 5/3; the relativistic correction is 
significantly larger in the case of the (p — n) model than in the (p — p) model. 
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